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I. INTRODUCTION 



The relativistic Toda lattice (RT) equation, 

qn =(1 + -qn-l){l + -gnj— ^ 



1 + ^exp{qn-i - qn) 



-{1 + -qn){l + -qn+i)——i : r, (1-1) 

c c 1 + ^exp(gn - gn+i) 

where qn is the coordinates of n-th lattice point, and " means the differentiation with 
respect to time t and c is the light speed, was introduced and studied by Ruijsenaars.^ 
This equation is derived by requiring the Poincare invariance to the well-known Toda 
lattice (TL) equation. The action-angle transformation and complete integrability were 
also investigated for the RT equation. Bruschi and Ragnisco^ constructed the recursion 
operator and Backlund transformation for eq. (1.1). The relation between the discrete 
time Toda lattice (DTL) and RT was found by Suris^ in the context of the Hamiltonian 
dynamics. We will show later that the DTL equation actually takes a part of the RT 
equation. In the limit of c — > cxo, the RT eq. (1.1) reduces to the ordinary TL equation, 

qn = exp(Q^_i - qn) - exp(Q^ - qn+i), (1.2) 

which is one of the most typical equations in soliton theory. 

The bilinear formalism introduced by Hirota^ is a key to understand the algebraic 
structure of soliton equations. By variable transformations, soliton equations are trans- 
formed into the bilinear form. The solutions for the bilinear equations are given in terms of 
certain determinants or Pfaffians. We analyze the RT equation from the view point of the 
bilinear theory. The Wronski determinant representation of solutions was first given for the 
KdV equation.^ Freeman and Nimmo^ have developed the Laplace expansion technique 
and given the Wronskian solutions for various soliton equations. Hirota, Ito and Kako^ 
showed that for the discrete analogues of soliton equations, solutions for their bilinear form 
are given by the Casorati determinant instead of Wronskian. 

M. and Y. Sato^'^ found that the soliton equations are reduced to the Pliicker relation 
with respect to the r function which is expressed by the infinite dimensional Wronski de- 
terminant. Date, Jimbo, Kashiwara and Miwa^°'^^ developed the theory of transformation 
group for solutions of soliton equations. Discrete soliton equations were also treated in 
their context. ■^■'^'■^^ Ueno and Takasaki-^^ extended Sato's idea to the Toda lattice hierarchy. 
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Recently Hietarinta and one of the authors (J.S.) found that the RT equation is 
transformed into a trihnear form through a suitable dependent variable transformation.^^ 
It was also shown that solutions for the trilinear equation are given by the determinant with 
Wronskian structure in two directions, rowwise and columnwise. The trilinear equation 
can be regarded as an extention of the bilinear form. Therefore it may be expected that 
the trilinear form has abundant algebraic structure. However the trilinear equation for 
RT is in the case that the coupling constant is negative, and hence the situation is not 
considered to be physical. To overcome this difficulty, we have investigated the solutions 
of eq. (1.1) and it has become clear that this equation holds various algebraic structures 
simultaneously. 

In this paper, we propose another type of determinant solution to the RT eq. (1.1), 
the Casorati determinant. This determinant represents the A'"-soliton solution for RT, so 
fits for the physical situation. The solution converges to the Wronskian solution for the 
one-dimensional TL equation in the limit of c ^ oo. The RT equation is transformed into 
three bilinear equations, the TL equation, the Backlund transformation (BT) for TL and 
the DTL equation. We show how the RT equation is composed of these three Toda systems. 
It is proved by using the Laplace expansion and the reduction technique that the Casorati 
determinant satisfies the bilinear equations. For discrete soliton equations, the difference 
formula is simply obtained by using the Casoratian technique. ""^^ Through the differential 
formula and difference formula for the Casorati determinant, the bilinear equations of RT 
are reduced to the "algebraic" identities. In II, the bilinearization of the RT equation is 
carried out by introducing the auxiliary dependent variables. In 1 1 1, we show the Casorati 
determinant solution for the bilinear equations which represents the A'"-soliton solution. A 
constructive proof that the determinant satisfies the RT equation is given in IV, and V. 
At first in IV, the solutions for the two-dimensional Toda lattice (2DTL) equation, BT 
of TL and discrete time two-dimensional Toda lattice (D2DTL) equation are derived in 
the form of the Wronski and Casorati determinants. Next in V, we explain the reduction 
technique which reduces the 2DTL and D2DTL equations into the TL and DTL equations, 
respectively. In VI, concluding remarks are given. 



II. BILINEAR FORM FOR THE RT EQUATION 
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The RT eq. (1.1) is transformed into the three bihnear equations, 

Dlfn-fn = 2{gngn-fn), (2.1a) 

{aDj; - l)fn ■ fn-l + 9n-l9n = 0, (2.1b) 

9n-l9n+l - fn^ O^Un+lfn-l - fl), (2.1c) 



through the variable transformations, 



gn=log%^, (2.2a) 

Jn 

t = X, (2.2b) 



c=- . (2.2c) 

In eqs. (2.1), g and g are the auxihary variables and D is the bilinear differential operator 
defined by^ 

D'j-g=id,-dy)'^fix)giy)\y=,. (2.3) 

As will be shown later, eqs. (2.1a), (2.1b) and (2.1c) are nothing but the bilinear forms of 
TL equation, BT of TL and DTL equation, respectively. 

Let us show that the RT equation actually follows from the bilinear eqs. (2.1). Dividing 
eqs. (2.1a), (2.1b) and (2.1c) by 2/^, fn-ifn and fn-ifn+i respectively, we get 

f^ = l + (log/„),„ (2.4a) 

Jnjn 

— = l + a(log-— )j„ (2.4b) 

Jn—ljn Jn 

9n-l9n+l 1,1-, 2\ fnfn (n a \ 

= a+{l-a)- , (2.4c) 

Jn—ljn+l Jn—lJn+1 



(2.5) 



where the subscript x indicates differentiation. We have an identity 

9n-l9n+l 9n9n _ 9n-l9n 9n9n+l 
fn—lfn+1 fnfn fn—lfn fnfn+1 

which is rewritten by using eqs. (2.4) as 

(a^ + (1 - a^)— + (log/„),,) = (1 + a(log ^).)(1 + a(log -Z^),). (2.6) 

In—lJn+l Jn Jn+1 

Noticing 

1 - o 



dx = " ^ ' " dt = —dt, 
c ac 



we get from eq. (2.6) 



1 11 
+ ^(log/n)« = (1 + -gn)(l + -gn+i)/(l + c^exp(gn+i - g^))- 

Finally subtracting eq. (2.9) from eq. (2.9) with n replaced by n — 1, we obtain 



(2.8) 



(2.9) 



-^Qn = (1 + -^n-l)(l+ -^n)/(l + C^exp(q'n - Qn-l)) 

- (1 + ^gn)(l + ^gn+i)/(l + c^exp(q'n+i - gn)), 



(2.10) 



which is the RT eq. (1.1) itself. Therefore if and Qn satisfy eqs. (2.1), defined by 

eq. (2.2a) gives the solution for RT equation. 

III. CASORATI DETERMINANT SOLUTION FOR THE RT EQUA- 
TION 



The Casorati determinant representation of soliton solution for the bilinear RT eq. 
(2.1) is given by 



and 



fn — 



Qn = 



(n 

(n 
^2 



(n 
4' 



in 



9n = 



(n+l) 



7(n+l) 



An) 



,(n) _ „-l 



Pi 



(n+N-l) 

(n+N-l) 
^2 



^(n+N-1) 
^(n+N-l) 



^(n+iV-1) 
^(n+N-l) 



1 -pitt 

Pi 



)"e^^ + ( 



P^ 



l-Pi^a 

Pi' 



-re' 



1 -Pia 



Te'^'+Pii )"e€' 



(3.1a) 



(3.1b) 



(3.1c) 



Pi a 



(3.2a) 
(3.2b) 



'^'^^^-P^^T^^^^'' +Pi\v^^r^^'^ (3-2c) 

?7i = piX + r]io, (3.3a) 

= Pi^x + ^io, (3.3b) 

where N is the number of sohtons, and Pi and rjio — cire the wave number and phase 
parameter of z-th soliton respectively. 

In the non-relativistic hmit, c tends to infinity, a becomes and 

9n+l = 9n-l = /n, (3.4) 

x^t. (3.5) 
In this case, eq. (2.1a) becomes the bihnear form of TL equation, 

DlU-fn = 2{fn+lfn-l-fl), (3.6) 

and eqs. (2.1b) and (2.1c) are trivial. Equation (3.6) is transformed into the TL eq. (1.2) 
through the variable transformation (2.2a). At c — > oo, the determinant fn in eq. (3.1a) 
gives the A^'-soliton solution to the TL equation in Wronskian form. 

We will prove by using the Wronskian technique and Casoratian technique that the 
bilinear eqs. (2.1) are satisfied by /„, and in eqs. (3.1). At first we give the determi- 
nant solution for the 2DTL, BT of TL and D2DTL equations in IV. Next in V, applying 
the reduction technique we recover the Casorati determinant solution in eqs. (3.1). 

IV. THREE TODA SYSTEMS AND THEIR DETERMINANT SOLU- 
TIONS 

In this section we show that the solutions for the 2DTL equation, BT of TL and 
D2DTL equation are given in terms of the Wronskian, Casorati determinant and also 
Casorati determinant, respectively. 

A. 2DTL 

The bilinear form of 2DTL equation is written as 

DxiDx_^Tn ■ Tn = 2(TnT„ - Tn+lTn-l). (4.1) 
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The Wronskian solution for this equation is given by^'-^^ 



(n) (n+1) (n+AT-l) 
(n) (n+1) (n+AT-l) 



(n) (n+1) (n+AT-l) 



(4.2) 



where (f^^^ 's are arbitrary functions satisfying the dispersion relations, 

din) (n+v) -, -, 

After Freeman and Nimmo,^ we use the following notation for simplicity, 



(4.3) 



|ni,n2, • • • ,nAr| 



(ni) (n2) 



(4.4) 



In this notation, in eq. (4.2) is rewritten as 



Tn = \n,n+ 1, - ■ ■ ,n + N - 1\ 



(4.5) 



We show that the above Wronskian satisfies the bilinear eq. (4.1) by using the Wron- 
skian technique. From eq. (4.3), the differentiation of the determinant |ni, n2, • • • , unI is 
given by 

N 

dxJni,n2,---,nN\ = ^\ni,n2,---,nj + iy,---,nN\, iy=l,-l. (4.6) 

i=i 



Thus we get the following differential formula for r^, 

dxiTn = \n,n + 1, ■ ■ ■ ,n + N - 2,n + N\, 
dx-i'Tn = |n — 1, n + 1, • • • , n + A'" — 2, n + A'" — 1|, 
dxidx_^Tn = \n-l,n-\-l,- ■ ■ ,n-\- N -2,n-\- N\ 
+ |n,n+l,---,n + A-2,n + A-l|. 

Now we consider an identity for 2 A x 2 A determinant, 

n-l I n+1 •■■ n+Ar-2 n+N \ n \ | n+AT-l 

n-1 10 1 n i n+1 " " " n+iV-2 1 n+AT-l 



(4.7a) 
(4.7b) 

(4.7c) 



= 0, (4.8) 
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where a number denotes a column vector of height N and means an empty block. 
Applying the Laplace expansion for determinant to the left-hand side of above equation, 
we obtain the algebraic bilinear identity for determinants, 



\n- l,n+l, - ■ ■ ,n + N -2,n + N\ \n,n+l, - ■ ■ ,n + N -2,n + N - 1\ 
-\n,n+l, - ■■ ,n + N - 2,n + N\ \n-l,n + l, - ■ - ,71 + N - 2,n + N -1\ 
+\n + 1, - ■ ■ ,n + N - 2,n + N - l,n + N\ \n - l,n,n+ 1, - ■ ■ ,n + N - 2\ ^ 0. (4.9) 

By using eqs. (4.7), this equation is rewritten into the differential bilinear equation for r^. 



(4.10) 



which recovers eq. (4.1). Hence we have proved that the Wronskian t„ gives the solution 
for 2DTL equation. 

B. BT of TL 

We have the bilinear form for BT of TL equation, -"^^ 



{aDx - 1)t'„ ■ Tn-l + TnT^_i = 0, 



(4.11) 



where a is a constant and x corresponds to xi in eq. (4.1), which transforms a TL Tn into 
another TL r^. Applying the BT recursively, we get an infinite sequence of TL's, r^, r^. 



T 



We denote the A;-th TL as t„(A;). Then the BT of TL is written as 

{aDx - l)Tn{k + 1) • Tn-l{k) + Tn{k)Tn-l{k + 1) = 0. 



The Casorati determinant solution for eq. (4.12) is given by 



(4.12) 



rn{k) = 



^^-\k) ^r'Hk) ... ^r''-'\k) 



(4.13) 



where <^^"^(A;)'s are arbitrary functions satisfying the dispersion relations. 



dx^\''\k) = ^^-^'\k), 



(4.14a) 



Ak4''\k) = 4''^'\k), (4.14b) 
where Ak is the backward difference operator with the difference intervai a defined by 

f{k)-f{k-l) 



For simplicity we introduce a convenient notation, 



(4.15) 



V^2''\kN) 



^^;:^\k,) ^r^(/c2) ••• <pT\kN) 



(4.16) 



In this notation, the solution for BT of TL, t„(A;) in eq. (4.13), is rewritten as 



Tn{k) = \nk,n+lk, ■ • • ,n + N - lk\, 



(4.17) 



or suppressing the index k, 



Tn{k) = |n, n + 1, • • • ,n + iV - 1|. 



(4.18) 



We show that the above Tn{k) actually satisfies the bilinear eq. (4.12). At first we 
investigate the difference formula for r. From eq. (4.14b) we have 



a<^J^+^)(/c + l) = <^f)(/c+l)-<^f)(/c). 

Noticing this relation, we get 

Tn{k + 1) = \nk+i,n+ lfc+i,n + 2fc+i, ■■■,n + N - 

= \nk,n+ lk+i,n + 2k+i, ■■ ■ ,n + N - lk+i\ 

where we have subtracted the 2nd column multiplied by a from the 1st column, 

= |nfc,n + lfc,n + 2k+i, ■ ■■ ,n + N - lk+i\ 

where we have subtracted the 3rd column multiplied by a from the 2nd column. 



(4.19) 



that is. 



= \nk,n+ Ik, - • • ,n' + N - 2k,n + N - lk+i\, 

Tn{k + 1) = |nfe,n+ Ife, • • • ,n + iV - 2/s,n + TV - lfe+i|, 
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(4.20) 



or suppressing the index A;, 



Tn{k+1) = \n,n+l,---,n + N -2,n + N - lk+i\. (4.21) 

Moreover in eq. (4.20), multiplying the A'"-th column by a and adding the {N — l)-th 
column to the A'"-th column, we get 

aTn{k + 1) = \nk, n+lk,---,n + N-2k,n + N- 2k+i\. (4.22) 

The important point is that the effect of difference appears only one column as shown 
in eqs. (4.20) and (4.22). Consequently the Laplace expansion becomes applicable. This 
technique for condensing the shift of index into only one column may be called a Casoratian 
technique. 

Differentiating eq. (4.22) with x and using eq. (4.14a) we obtain 

adxTn{k + 1) = \nk,n+ Ik, ■ ■ ■ ,n + N - 3k,n + N - lk,n + N - 2k+i\ 
+ \nk,n+ Ik, - ■ ■ ,n + N - 3k,n + N - 2k,n + N - lk+i\ 
= \nk, n + Ik, - ■ ■ ,n + N - 3k,n + N - lk,n + N - 2k+i\ 
+ Tn{k+1), (4.23) 

where we have used eq. (4.20) in the last line. We have also 

dxTn-i{k) ^\n-l,n,---,n + N -3,n + N (4.24) 

Let us introduce an identity for 2N x 2 A'' determinant, 

Tifc ■■■ n+N-3k I n+N-lk \ n+N-2k+i \ ® \ n+N-2k 

i n+JV-lfe I n+JV-2fc+i I n-lfc nu ' " n+N-3k '' n+N-2k 



= 0. (4.25) 



Applying the Laplace expansion to the left-hand side, we obtain the algebraic bilinear 
identity for determinants, 

Infe 5 •■■ ) n+N—3k , n+N—l^ i n+Af — 2te_|_i | |n— Ife ; n^, n+N—3k i n+AT— 2^ | 
~\nki ■■■■> n+N—3k i n+N—2k ) n+JV— 2fe+i | In— Ifc ,nkJ ■■•■> n+N—3k ) n+N—lk I 
— Ufc, n+AT-Sfe, n+AT-lfc, n+JV-2fe lln-lfe, nfc, n+JV-3fe5 n-|-Ar-2fc+i | = 0, (4.26) 
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which is rewritten by using eqs. (4.20), (4.22), (4.23) and (4.24), into the differential bihnear 
equation, 

(ada^Tnik + 1) - Tnik + l))Tn-iik) - ar^ (/c + l)aa;Tn_i (/c) + t„(A;)t„_i (/c + 1) = 0. (4.27) 

This equation is equal to eq. (4.12). Hence we have proved that the Casorati determinant 
Tn{k) in eq. (4.13) gives the solution for the BT of TL equation. 



C. D2DTL 



Using the Casoratian technique, we show that the D2DTL equation is satisfied by the 
Casorati determinant. The bilinear form of D2DTL equation is given by^'^^ 



Tn{k + lJ + l)Tn(/c, - Tn{k + 1, l)Tn{k, I + 1) 
--ab{Tn{k + 1,1 + l)Tn{k, I) - Tn+l{k + 1, l)Tn-l{k, I + 1)), 



(4.28) 



where depends on two discrete independent variables k and /, and a and b are the 
difference intervals for k and / respectively. 

The Casorati determinant solution for eq. (4.28) is given as 



Tn{k,l) = 



^^;:\k,i) ^Sv"+^^(/c,o 



(4.29) 



where <^^"^ 's are arbitrary functions of k and I which satisfy the dispersion relations. 



Akvi^\k,l) = 4''+'\k,l), 

Here and Ai are the backward difference operators defined by 

f{k,l)-f{k-l,l) 



(4.30a) 



(4.30b) 



AkfikJ) 



Aif{k,l) = 



a 



f(k,l)-f{kj-l) 



(4.31a) 
(4.31b) 
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We use the following simple notation similar to that in eq. (4.16), 



<PN''\kN, In) 



(4.32) 



In this notation, the solution to D2DTL equation in eq. (4.29) is 



Tn{k,l) = |n^,n+l^,---,n + iV-l^ 
I I I 

Neglecting the index k and I, we write 

Tn{k,l) = \n,n + 1, ■ ■ ■ ,n + N — 1\. 



(4.33) 



(4.34) 



Now let us examine the difference formula for the above Tn- Note that the dispersion 
relation for k in eq. (4.31a) is equal to the one for the solution to BT of TL in B. Hence 
in the same manner as shown in eqs. (4.20) and (4.22), we obtain 



Tn(/c + 1, /) = |n ^ , n + 1 ^ , • • • , n + iV - 2 ^ , n + iV - 1 
/ I I 

aTn{k + l,l)= \n^,n + l^,---,n + N -2^,n + N -2 
Similarly we get 



/ / I I 



(4.35) 
(4.36) 



Tn{k,l + l)^\n ^ ,n + l^,---,n + Ar-2^,n + iV-lJ, 



bTn{k,l + l) = \n + l ^ ,n+l^,---,n + N -2^,n + N -1^ 

l+i I II 



(4.37) 
(4.38) 



because the role of I and b is parallel to that of k and a except that the ordering of index 
n is reversed. In short, we write 



Tnik+l,l)= \n,n+l,---,n + N -2,n + N -lk+i\, 

aTn{k + 1,1) = \n, n + 1, ■ ■ ■ , n + N — 2, n + N — 2k+i\, 
Tn{k,l + 1) = \ni+i,n + l,---,n + N -2,n + N 
bTn{k,l + l) = \n+li+i,n + l,---,n + N -2,n + N -1\ 
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(4.39) 

(4.40) 
(4.41) 
(4.42) 



We give the Casoratian technique for two variable shifted case. From eq. (4.41), 
Tn{k + l,l + 1) is given by 

Tn{k + l,l + l) = \n^_^^,n+ Ife+i, ■■■,n + N- 2k+i,n + N- lk+i\. (4.43) 

In the following we show that the shifts of index k is condensed into only the most right 
column of the determinant. Prom eqs. (4.30), v?f ^ satisfies 

AkA^i^^ = (4.44) 

that is, 

^^\kj) - ^t\k,l- 1) - ^t\k - 1,0 + ^t\k - 1,/ - 1) = ah^t\k,l), (4.45) 
which is rewritten as 

(1 - ah)^f {k + 1,1 + 1) = (/c + 1, + {k,l + l)- {k, I) 

= )(/c, l + l) + a^i^+'\k + 1, 0- (4.46) 

Therefore multiplying the both hand sides of eq. (4.43) by (1 — ab) and rewriting the 1st 
column of the determinant by the use of eq. (4.46), we obtain 

{l-ab)Tn{k + 1,1 + 1) 

= |n ^ ,n+lk+i,---,n + N -2k+i,n + N -lk+i\ 
l+i 

+ a|n + 1 , n + U+i, • • • , n + TV - 2fe+i, n + N- lk+i\. (4.47) 
I 

The second term in r.h.s. vanishes because its 1st and 2nd columns are the same. So we 
get 

(1 - ab)Tnik +1,1 + 1) = \ni+i,n + l^+i, • • • , n + TV - 2k+i,n + N- lk+i\. (4.48) 

In the above determinant, by subtracting the (j + l)-th column multiplied by a from the 
j-th. column for j = 2, 3, • • • , A?" — 1, Tn{k + 1, Z + 1) is given by 

(1 - ab)Tn{k+ 1,1 + 1) = \ni+i,n+l,---,n + N -2,n + N - lk+i\. (4.49) 
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As is shown, even if the two variables k and I are shifted, t„ is also expressed in the 
determinant form whose columns are almost unchanged and only edges are varied. Thus 
we can use the Laplace expansion technique. 

We consider an identity for 2N x 2N determinant, 





nj+i I n+1 



n+N-2 n+N-lk+i ' n 







n I n+1 



n+N-1 
n+N-2 ! n+N-1 



= 0. 



By the Laplace expansion, we have 



(4.50) 



|n;+i ) n+1) •••) n+N-2 J n+N-l^+i I In; n+1) ■••) n+N-2j n+N-1 \ 

~\m n+1) •••) n+Ar-2) n+iV-lfe + i Hnj + D n+1) •••) n+N-2i n+Af-l| 

+ |n+l) •••) n+N-2, n+N-1, n+N-l^+i I |n;+i) n) n+1) ■••) n+JV-2| = 0. (4-51) 

Substituting eqs. (4.39)-(4.42) and (4.49) into eq. (4.51), we obtain 

(1 - ab)Tn{k +1,1+ l)Tn{k, I) - Tn{k + 1, l)Tn{k, I + 1) 

+ aTn+i{k + 1, l)bTn-i{k, Z + 1) = 0, (4.52) 

which recovers eq. (4.28). This completes the proof that the Casorati determinant is indeed 
the solution for D2DTL equation. 

V. REDUCTION TO THE RT EQUATION 



Summarizing the arguments in IV , we see that the three bilinear equations, 

DxDyTnik, I) ■ Tn{k, I) = 2(Tn(/c, /)Tn(/c, /) - Tn+i(/c, /)Tn_i(/c, /)), (5.1a) 
(aDa, - l)Tnik + 1,1) ■ Tn-lik,l) + TnikJ)Tn-lik + 1,1) = 0, (5.1b) 



Tn{k +1,1+ l)Tn{k, I) - Tn{k + 1, l)Tn{k, I + 1) 
= ab{Tn{k +1,1 + l)Tn(k, I) - Tn+l{k + 1, l)Tn-l{k, I + 1)), 



are simultaneously satisfied by the Casorati determinant. 



Tn{k,l) 



^^-\k,l) ^t^'\k,l) 
^^^\k,l) ^^^+'\k,l) 



^r^'-'HKi) 

(n+N-1) 
V2 



{k,l) 



v^;:\k,i) v^^^'\k,i) ... ^S^+^-^^(/c,o 
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(5.1c) 



(5.2) 



supposed that ^'i'^^'s obey the dispersion relations, 



o (n) (n+l) 


(5.3a) 




(5.3b) 




(5.3c) 




(5.3d) 



where and Ai are the same as those in eqs. (4.31). For example, we can take (/p^"^'' as 
=pr(l - -h)-'e^^+q^{l - q,a)-\l - -6)-'e^^ (5.4) 

Pi Qi 
1 

r]i = —y + piX + Tjio, (5.5a) 
Pi 

Ci = —y + QiX + ^io, (5.5b) 

Qi 

where pi, qi and rjio, ^io are arbitrary constants which correspond to the wave numbers 
and phase parameters of solitons, respectively. It can easily be seen that cp^"'^ in eq. (5.4) 
actually satisfies eqs. (5.3). 

The reduction is the technique to derive subhierarchies of equations by restricting the 
solution space, in other words, imposing some conditions to the solutions. In the case of 
RT equation, we take in eq. (5.4) 

Qi = —, (5.6) 
Pi 

and 

b = a. (5.7) 

On the conditions (5.6) and (5.7), ip\"'^ in eq. (5.4) satisfies 

(a, + a>f) = (p, + l)v.f), (5.8a) 

Pi 



^f\k - 1,Z - 1) = {l-Pia){l - -a)^t\k,l). (5.8b) 

Pi 



Hence we obtain 



{dx + dy)Tn = Stu, (5.9a) 

Tn{k-l,l-l) = PTn{k,l), (5.9b) 
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where 



^ 1 

P = l[{l-pia){l--a). 



1=1 



By using eqs. (5.9), eqs. (5.1a) and (5.1c) are rewritten as 

DlTn{k, I) ■ Tn{k, I) = 2{Tn+l{k, l)Tn-l{k, I) - Tn{k, l)Tn{k, I)), 



(5.10a) 
(5.10b) 

(5.11a) 



Tnik, l)Tn{k, I) - Tn{k + 1, l)Tn{k - 1, I) 
= 0?{Tn{k, l)Tn{k, I) - Tn+l{k + 1, l)Tn-l{k - 1, I)). 

Here we may drop y and /-dependence. Thus eqs. (5.11) and (5.1b) give 

DlTn{k) ■ Tn{k) = 2{Tn+l{k)Tn-l{k) - Tn{k)Tn{k)) , 

(aDx - l)Tn{k) ■ Tn-l{k - 1) + Tn{k - l)Tn-l{k) = 0, 
Tn{k + l)Tn{k - 1) - Tnik)Tn{k) = tt^ {Tn+l{k + l)Tn-l{k - 1) - Tn(/c)Tn(/c)). 



(5.11b) 



(5.12a) 

(5.12b) 
(5.12c) 



Equations (5.12a) and (5.12c) are nothing but the bihnear form of the one-dimensional TL 
and DTL equations. From eqs. (5.2) and (5.4)- (5.7), the Casorati determinant solution of 
eqs. (5.12) is 

^[-\k) ^r'\k) ... ^r^'-'Hk) 
^^-\k) ^r'\k) ... ^r^-^^(/c) 



Tnik) 



Letting 



Pi Pi 

r]i =piX + r]io, 

Pi 



fn = Tn(n), 
9n = Tn-l(n), 
9n = Tn+l(n), 
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(5.13) 

(5.14) 
(5.15a) 
(5.15b) 

(5.16a) 
(5.16b) 
(5.16c) 



eqs. (5.12) with k = n become the bihnear form of RT (2.1), and eq. (5.13) gives its solution 
(3.1). Thus we have completed the proof that the A'"-soliton solution for the RT equation is 
given by the Casorati determinant in eqs. (3.1). We note that besides the soliton solution, 
other solutions such as the rational exponential soliton can be obtained in the same way. 



VI. CONCLUDING REMARKS 



We have shown that the solution of the RT equation is given by the Casorati determi- 
nant. The composition of three Toda systems derives the relativistic modification of the 
TL equation. This fact means that the Toda equation has quite fruitful content. 

The difference between the relativistic and non- relativistic TL's is only the choice 
of dependent variables in the D2DTL system. In the two-dimensional space of discrete 
independent variables n and k, we have the Casorati determinant solution Tn{k). If we 
take rn{n) as the dependent variable we obtain the relativistic version of TL through 
the variable transformation (2.2a). The difference interval a gives the light speed c. On 
the other hand, adopting Tn{k) with fixed k as the dependent variable, we get the original 
TL in the same variable transformation. Namely, the RT is a sublattice of the D2DTL 
while the TL is another sublattice. 

It is well known that the other soliton equations such as the KdV, modified KdV, 
sine-Gordon and nonlinear Schrodinger equations are closely related to or derived from the 
Toda equation by suitable limiting procedures. Therefore we expect that the relativistic 
version of these equations can also be derived from the RT equation. We conjecture how 
to construct such relativistic soliton equations and their solutions. 

(1) For a given soliton equation, consider its integrable discretization. 

(2) Derive the proper two-dimensional lattice system (before reduction). 

(3) Take the different sublattice, that is, the different direction of propagation. 

(4) In a suitable variable transformation, the relativistic equation corresponding to the 
original one will be obtained. 

It may be interesting to investigate the relativistic version of nonlinear integrable equations 
in this scenario. 

It should be noted that the RT equation has many other aspects and there are several 
types of determinant solutions besides the Casorati one. Actually the Wronskian of two- 
directional type yields quite different solutions from the present A/"-soliton solution. We can 
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derive some interesting solutions from the Wronskian, which will be reported elsewhere. 
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